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The pre-buckling part of this analysis is based on Lekhnitskii’s (ref. 3) complex variable 
equations. For a two dimensional stress analysis, the equilibrium equations are: 



The solutions for the stresses in terms of a function F (Airy stress function) are: 
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Substituting the Laurent series representation into the force equation and evaluating from i to i-1 
gives: 

(12 & 13) 
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where the Akn are the unknown Laurent series constants, Cmkn contains the coefficients of the 
Laurent series constants, and F m are the resultant applied forces or applied displacements. To 
improve the solution of this system of equations, a least squares approach is taken. Using this 
method, twice as many equations will be used as there are unknowns. Therefore: 



{ Ak n } is a 8 *N vector 
{ Fm } is a 16 *N vector 

To solve this system, each side of the system is multiplied by [ Cmkn as follows: 
[ Cmkn [ Cmkn ] { Akn } = [ Cmkn 1^ { Fm } 
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Net-Section Shear Stress for a Square Plate with 5% Hole 
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Comparison of Net-Section Shear Stresses with Finite Element Results 
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